Mesoscopic Aspects of Strongly Interacting Cold Atoms 
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Harmonically trapped lattice bosons with strong repulsive interactions exhibit a superfluid-Mott- 
insulator heterostructure in the form of a "wedding cake". We discuss the mesoscopic aspects of 
such a system within a one-dimensional scattering matrix approach and calculate the scattering 
properties of quasi-particles at a superfluid-Mott-insulator interface as an elementary building block 
to describe transport phenomena across such a boundary. We apply the formalism to determine the 
heat conductivity through a Mott layer, a quantity relevant to describe thermalization processes 
in the optical lattice setup. We identify a critical hopping below which the heat conductivity is 
strongly suppressed. 

PACS numbers: 05.30.Jp, 74.78.Na 



I. INTRODUCTION 

Cold bosonic gases subject to an optical lattice allow 
for an accurate emulation of the Bose-Hubbard model 
of interacting lattice bosons,^ with a short range (on- 
site) interaction between bosons, a hopping restricted to 
nearest- neighbors, and no inter-band transitions, at least 
for sufficiently deep lattice potentials. Furthermore, the 
system is almost perfectly decoupled from environmen- 
tal degrees of freedom during typical experimental times. 
These aspects make cold atoms a perfect test bed for 
probing lattice Hamiltonians relevant to condensed mat- 
ter physics .- 

There is one feature, however, where cold atom systems 
differ from the generic solid-state setup, as they typically 
experience an inhomogeneous potential due to the trap, 
and thus inferring bulk properties is often hampered by 
finite size effects<^i^ On the other hand, as a result of this 
confinement, new interesting structures and effects may 
occur, with the "wedding cake" involving layers of Mott- 
insulating and supcrfluid phases of strongly correlated 
bosons providing a prominent example^i^ [cf. Fig. [TJa)]. 
The evolution of the ground state (gs) across this inho- 
mogeneous system can be easily understood within the 
framework of a local density approximation combined 
with a mapping between the position in the trap and 
the corresponding point in the bulk phase diagram;^ The 
most pronounced new feature in this layered structure 
is the superfiuid-Mott-insulator (S-MI) interface. Such 
boundaries between phases with different symmetries are 
known to exhibit interesting effects; a well known exam- 
ple is the phenomenon of Andreev reflection^, between a 
normal-metal and a superconductor, where electrons in- 
cident on the superconductor from the normal metal are 
reflected back in the form of a hole retracing the elec- 
tron's path. Given the strongly correlated nature of the 
two phases framing the S-MI interface, similar interesting 
phenomena may be expected in the present case. In this 
work, we provide a description of the scattering proper- 
ties of elementary excitations at the S-MI interface. Such 
knowledge then allows us to determine the energy (or 



FIG. 1: (a) One-dimensional model of a superfluid-Mott- 
insulator (S-MI) heterostructure in a trap defined by the po- 
tential Vi at position i. The dark areas denote Mott-insulating 
regions with a fixed density pi, while the bright areas are su- 
perfluids of condensed particles (center, Sp) and condensed 
holes (wings, S^). (b) Sketch of a scattering event at the 
superfluid-Mott-insulator interface, where a phonon incident 
from the superfluid generates back-reflected quasi-particles of 
phonon- and amplitude-type as well as particle- and hole-type 
excitations transmitted into the Mott insulator. 



heat) flow across interfaces and we determine the heat 
conductivity across the Mott insulator region connecting 
two superfluids. 

Mesoscopic aspects in a wedding cake structure have 
been studied by Vishveshwara and Lannert^ who cal- 
culated the Josephson coupling across a Mott-insulator 
domain (a S-MI-S junction) supported by the over- 
lap between exponentially suppressed ground state wave 
functions (superfluid order parameters) across the Mott 
phase. Here, we are interested in the scattering dy- 
namics of the excitations near a S-MI boundary and 
the transport associated with them across one or mul- 
tiple interfaces; these excitations are sound (Goldstone) 
and massive (Higgs) modes within the superfluid phase, 
and particle- and hole-type excitations in the Mott 
insulatorii^ 

After developing the general framework describing 
transport in an inhomogeneous system with phase 



boundaries, we calculate analytically the reflection and 
transmission coefficients for a phonon mode of the super- 
fluid incident on a Mott insulator [cf. Fig. [Hb)] . While 
these coefficients can be used to describe the scattering 
of a wave-packet of phonons (a density disturbance), here 
they mainly serve as an illustration of our approach. We 
then use our results (including those for scattering of 
massive modes) to calculate the heat transport through 
a Mott barrier. Knowledge of the heat conductivity k al- 
lows to estimate the thermal contact between superfluid 
shells in the wedding cake structure of a trapped Bose 
system. Note that temperature (actually entropy) im- 
balances quite naturally occur in optical lattice systems 
as the lattice potential is turned on and entropy is ex- 
pelled from the newly formed Mott-insulating regions.-- 
Our heat conductivity k then relates to the thermaliza- 
tion process across different superfluid rings. 

Before developing our formalism in detail, we give an 
overview of the ideas and concepts utilized in this pa- 
per. The excitations close to one of the Mott-insulating 
lobes involve coherent superpositions of various site oc- 
cupation numbers, and their wave functions are given 
by a four-spinor structure. This is reminiscent of the 
two-spinor structure of excitations in the Bogoliubov- 
de Gennes equations describing an inhomogeneous su- 
perconductor (note that for most unconventional super- 
conductors, the reduction to a two-spinor is possiblc^^). 
The program to be carried out in order to flnd the trans- 
mission, reflection, and transformation of quasi-particles 
then is identical to the one introduced by Blonder, et 
alJ^ for the superconductor-normal-metal boundary (in 
order to simplify the analysis, we consider here a one- 
dimensional situation and leave geometric effects due to 
finite impact angles for a later study): First, we deter- 
mine the excitation energies tvjJxik) and the four-spinor 
structure X(A;) for the quasi-particle excitations in the 
superfluid (X = s, m: sound and massive modes) and in 
the Mott-insulator (X = p, h: particle and hole modes); 
below, this will be done for a translation invariant sit- 
uation in a second-quantized formalism. Second, we 
switch to a first-quantized formulation and account for 
the inhomogeneous setup; for a slowly varying potential 
\Vi -^ V{y)] this can be done within a quasi-classicalii 
or Wentzel-Kramers-Brillouin (WKB) ^^i^^i^^ approxima- 
tion, with the wave functions assuming the form 



assume the form 



^x(y) = 



1 



V\k%M 



exp 



t / dxkj,{x) X[kj,{y)]. (1) 



For a quasi-particle with energy e, the wave vector 
kxiy) is obtained via proper inversion of the disper- 
sion huJx{k)[V{y)] = e, where the potential V{y) en- 
ters the expression via a shift of the chemical potential, 
S^ ^ 6^ — V{y). Third, we account for different phases 
in the setup by calculating the transfer matrix across the 
interface^- via matching of the wave functions and their 
derivatives at the boundary. For the case of a phonon 
incident from a (particle-type) superfluid Sp on a Mott- 
insulator MI [cf. Fig. mb)], the wave functions locally 






(2) 



with ip^ (y) denoting plane-wave functions of type ([T|) 
with constant wave vector ±fc. The scattering ampli- 
tudes Tpg (transmitted particle), r^^ (transmitted hole), 
rss (reflected sound), and rms (reflected massive mode) 
are obtained from the continuity conditions across the 
interface, t/>^(0) = ^^(0) and ^y^Pl{0) = dyxPlM- 

In order to fully characterize the Sp-MI boundary, 
one has to determine the 4x4 transfer matrix relat- 
ing phonon and massive modes propagating to the left 
and right through the superfluid Sp with the particle and 
hole modes propagating to the left and right through the 
Mott-insulator, requiring the solution of four scattering 
problems of the above type [Eq. ([2])] (the corresponding 
task has to be solved in order to describe the interface be- 
tween a hole- type superfluid Sh and a Mott insulator). In 
the following, we concentrate exclusively on the scatter- 
ing properties at the interface where the superfluid order 
parameter V' vanishes; this guarantees the solvability of 
the matching conditions at the boundary. As our final re- 
sult, we present: (i) simple analytical expressions for the 
scattering coefficients at small hopping i ^ and near 
the critical value tc at the tip of the Mott-insulator lobe, 
as well as numerical results for two cases in between; and 
(ii) the behavior of the heat conductivity k as a function 
of the model parameters and the temperature T. 

The body of the paper is organized as follows: In 
Sec. ini we derive the elementary excitations in the vicin- 
ity of the Mott-insulating regions of the Bose-Hubbard 
model. We briefiy review the mapping to a spin-1 
problem-'^*' before discussing the spinor structure and 
symmetry properties of the wave functions. Section [TTTl is 
devoted to the calculation of the scattering coefficients for 
a superfiuid-Mott-insulator interface, and the heat con- 
ductivity through a Mott barrier is determined in Sec. 
TV\ We summarize our results and conclude in Sec. IVl 



II. EXCITATIONS 

In this section, we first derive the excitations of the 
strongly correlated superfiuid and the Mott-insulating 
phase and exploit the time-reversal symmetry of the 
problem to obtain a well-suited set of spinor wave func- 
tions. For the reader who is less interested in the tech- 
nical details, the energies fvjjx{k) in Eqs. Q and (O, 
the spinor wave functions X(fc) in Eqs. ^ and ([7]), and 
their time-reversal invariant combinations [Eq. (J12p ] , rep- 
resent the main result of this section and Table |T] provides 
a physical interpretation of the four-spinor components. 

We find the dispersion and eigenfunctions of quasipar- 
ticle excitations of the Bose Hubbard- model following the 
procedure described in detail in Ref. ilCl Starting with 
the Bose-Hubbard Hamiltonian (the bosonic operators a^ 



create particles in Wannier states at site i, Siii = ajoi— no 
measures the density from a mean density uq, t is the 
hopping energy, U accounts for the on-site interaction, 
5^ is the chemical potential controlling the particle num- 
ber, and Vi accounts for the harmonic confinement), 



H 






we first consider a homogeneous situation with Vi = 
and truncate the bosonic Hilbert space to a site basis 
with three local states \0)i, |±l)i- The state \0)i refers to 
no bosons on site i, while the states [a)i include one more 
(less) particle. Within this restricted space, we assume 
a trial wave function for the ground state |^) — Yii \i^i), 
where 

|V',)=cos(7?)|0),+sin(^)[sin(x)| + l)., + cos(x)|-l).]. 

Minimizing the variational energy e„^,. ~ {'^111^^1'^) with 
respect to cr = 7r/4 — x and subsequently expanding e^^^ 
in the order parameter tp = (^|ai|^) = •\/rio/2sin(2^), 
we obtain e^„ « aip'^ + bip'^/2. To simplify expressions, 
we assume large filling with -yno"+T « ^/no- With 
the critical hopping tc = U/Suq, we find for the order 
parameter close to the upper (Sj,-MI) and lower (S/j- 
MI) phase boundaries <5/i^ — ±[/-y/l — t/tc/2 the result 



Af Jt, where A±^ = Sfif ± Sfi. 



The calculation of the excitations above the mean- 
field ground state |^) involves a spin- wave analysis in 
a slave boson description; the ground state in this ef- 
fective pseudo-spin 1 formalism defines the direction of 
the magnetic order, while the remaining two degrees 
of freedom provide the excitations. For i/) — > 0, the 
result can be expressed as a 4x4 matrix Hamilto- 
nian H = J2mn T^T^mnTn with four-spinor operators 
T = {tlj,,t^_^i^,ti^^k,t-i-k) describing particle and 
hole creation and destruction above a mean-filling uq 
(cf. Table H]) (the operators t]^^ ^ act on the vacuum 
states |vac)i and relate to the bosonic operators aj via 
ti,^, - (at)»"±VV(^^^±T)!). 

The diagonalization via a Bogoliubov transformation 
provides the dispersions (we measure wave vectors k in 
units of 1/a, with a the lattice constant) 

U 



f^Pih) (k) = — \/l - (t/tc) cos(fc) =F Sfi (4) 



in the Mott-insulator and 



^sM(fc) = 4inoy'a2^„)-/32^„) 



fc^o J ck, 
|A„ 



(5) 



in the superfiuid phase for -0 —*■ at i < ic, where 

as = [2/(2-i/ic)-cos(fc)]/2, a„ = 2i,/i -^ a„ 

Ps = -pm = tcos{k)/{2U-t). 



four-spinor 


Microscopic interpretation 


Dirac 




r i,, ] 


create particle 


e-,T 




^u. 


create hole 


e-,i 




ti,-k 


destroy particle fluctuation in gs 


e+,T 




U-i,-J 


destroy hole fluctuation in gs 


e+,i 



TABLE I: Physical interpretation of the four-spinor and its 
similarity to the Dirac spinor. In both cases, we deal with 
two species of "particles": the particle and hole excitations 
discussed here correspond to spin-up/down electrons in the 
Dirac equation. The role of anti-particles in the Dirac spinor 
is played by ground state fluctuations in the present context: 
the anti-particle of a particle is a fluctuation with a "missing" 
particle in the gs and similar for the hole. Hence the top and 
bottom entries are both particle-type (creation of particle and 
destruction of a hole), while the two middle entries are both 
hole-type (creation of a hole and destruction of a particle). 



The sound velocity c of the phonon is c = \/4Utno 
/{2tc/t—l) and the gap A^ of the massive (or amplitude) 
mode is given by A^ = C/^1 — t/tc- At the Sp-MI inter- 
face, the particle excitation in the Mott-insulator trans- 
forms into the (particle- type) sound mode in Sp, while 
the hole excitation in the Mott-insulator transforms to 
the (hole- type) massive mode in Sp. Correspondingly, the 
bottom of the hole branch in the Mott-insulator matches 
up with the gap A^ of the massive (hole-type) mode 
in Sp, while the bottom of the particle branch in the 
Mott-insulator goes to zero. The use of Eq. ([T]) within 
an inhomogeneous superfiuid phase requires generaliza- 
tion of the result in Eq. (O to finite values of tp, which 
provides the dependence on 6fi away from Sfj,^ and hence 
on the local potential V{y). Within the Mott region, the 
particle- and hole spectra undergo a simple shift ±6fi [ 
cf. Eq. ^], with a corresponding simple dependence on 
the smooth potential V{y). 

The spinor eigenstates in the Mott insulator are 




Pk = S n } and 




(6) 



with the coefficients 



A{k) = cosh(atanh{(t/tc)cos(fc)/[2-(f/ic)cos(A:)]}/2), 
B{k) = sinh(atanh{(i/ic)cos(fc)/[2-(f/ic)cos(A:)]}/2), 

providing the "dressing" of a particle (with amplitude A 
in pfc) by missing hole- type fluctuations (with amplitude 
B) in the ground state (and vice versa for h^) and fulfill 
the normalization condition. 

The spinor eigenstates in the particle-condensed super- 



fluid phase Sp are 

n^As{k) 
-^+Bs{k) 

n^Bs{k) J 

with the cocfBcients 



and m^ 



n+A,n{k) 
n-B^{k) 
n+B„^{k)} 



(7) 



^s(m) =cosh{ata-nh[Ps{m){k)/as(„i){k)]/2}, (8) 
Bs{„i) =sinh{atanh[/3^(„,)(fc)/as(„)(A:)]/2}. (9) 

Furthermore, Q± ^ [cos(/(i)) ± sin(/(t))]/\/2 with 
f{t) = arctan[2(ic/i)-\/l — t/tc]. The eigenstates for S/i 
are obtained by replacing 



n. 



n^. 



(10) 



The nature of these excitations is easily understood 
in the limits t — » and t ^ tc- For t — > 0, the coeffi- 
cients A, ^s(m) -^ 1 and B, i?s(m) ^ 0, telling us that the 
ground state turns into a classical one devoid of fluctua- 
tions [cf. Table I]. The excitations [Eq. ^] in the Mott 
phase combine particles with absent hole fluctuations and 
holes with absent particle fluctuations and hence are of 
particle and hole types, respectively. With f7+ -^ 1 and 
51- -^ the weights are set differently in the superfluid 
phase; here, the excitations involve particles dressed with 
holes and holes dressed with particles, reflecting the col- 
lective nature of these excitations. Crossing the bound- 
ary (5/1+ into the particle- type superfluid Sp, the particle 
mode of the Mott insulator condenses, imprinting the 
predominant particle nature onto the sound mode [the 
gapped hole mode goes over to the massive (amplitude) 
mode in the superfluid]. 



For t 



tr, the coefficients A, A 



B.B 



s{v\) 



and 



s(jn) 



oo; the excitations exploit the strong fluc- 
tuations in the ground state but keep their particle and 
hole character in the Mott insulator. With ^± — > l/\/2, 
both sound and massive modes have a mixed particle-hole 
character and draw large weights from the ground-state 
fluctuations. 

In the next section, we will have to match the spinor 
wave functions and their derivatives, providing us with 
eight conditions for only four unknown scattering coeffi- 
cients. In order to eliminate the additional spurious con- 
ditions, it is convenient to exploit the symmetry under 
time reversal and we introduce the time-reversal operator 



T = C 




(11) 



where C denotes the action of complex conjugation and 
ctq is the identity matrix in two dimensions. The form 
of T can be inferred from the behavior of a^ under 
time reversal and the relations between the operators ai 
and ii,Oj^i,±i- In addition to the above spinors X = 



p, h, s, m, we define the four linearly independent eigen- 
vectors TX. For the situation with unbroken T invari- 
ance discussed here^^ it is convenient to define the new 

states, 



^i 



(XfeTiTXfc)/V2. 



(12) 



These form a basis of T representations and are multi- 
plied with ±i under the operation T. Without loss of 
generality, we choose to work with the "-I-" -eigenstates 
and drop the superscript in the following. 

With these spinors at hand, we are now in the posi- 
tion to describe inhomogeneities. For piecewise constant 
"potentials" , we can use plane- wave-type left and right 
moving spinor wave functions 



Vf"(2/) = e*(^+'^)/4e±''=-«Xfe., 



(13) 



and combine them into the scattering states defined in 
Eq. ^. For slowly varying parameters ("potentials"), 
the states [Eq. p^ ] determine the spinor structure in 
Eq. H]). 



III. SCATTERING COEFFICIENTS 

Next, we discuss a specific situation, the scattering of 
a phonon on a Mott-insulator boundary, and determine 
the scattering coefficients by matching the wave func- 
tions [Eq. (21)] at the boundary between the superfluid 
and the insulating region. The scattering state ([2]) de- 
scribes a phonon excitation incident from a particle-type 
superfluid Sp onto a Mott insulator [cf. Fig. [1]. These 
excitations are relevant in an experiment where a den- 
sity perturbation is applied to the middle of the trapj^ 
while one has to describe the density perturbation by a 
wave packet of phonons, here, we only discuss the scat- 
tering properties of plane- wave excitations as the generic 
building block. Imposing the continuity conditions across 
the Sp-MI interface, we obtain the scattering coefficients 
^ss, rms, Tps, and Ths', only the first two spinor elements 
are relevant as the second two satisfy the matching condi- 
tions automatically due to T invariance. In the following, 
we discuss these coefficients in the limiting cases t — > 0, tc 
and illustrate their behavior at intermediate values of t 
in Fig. [21 Furthermore, here, we restrict the discussion 
to those modes propagating in the superfluid as well as 
in the Mott-insulator region; in the next section, where 
the heat transport through a finite Mott region is dis- 
cussed, the contribution of evanescent modes has to be 
considered as well. 

For i — > 0, a sound mode incident from Sp can be 
reflected as a sound mode or transmitted into the Mott 
insulator as a particle excitation; hence, the coefficients 
connecting to propagating scattered modes are r^s and 
Tps] given the energy e of the incident phonon mode, they 
can be conveniently expressed through the wave vectors 
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FIG. 2: Scattering coefficients for an incoming sound mode at an energy e = fiajs(7r)/4. (a) For all values of t/tc- The full 
lines show the exact values, while the dash-dotted line is the approximate value from Eq. (|14p . (b) Comparison between the 
exact and approximate values for ras- (c) All coefficients around the "critical" value t ~ 0.971 where for e = /iaJs(7r)/4 the 
hole and massive modes become exponentially damped. Note, that values of the scattering coefficients larger than 1 are not in 
contradiction with particle conservation as these excitations do not have a well defined particle character, (d) For values i « f^ 
the dash-dotted results from Eq. (|16fl again describe the exact solutions well. 



kl and k^ of the sound- and particle excitations involved, 
ki - K 






^1 + ^p 
k% -t- ki, 



t cos{kl) 



1- 



4ic 1 - cos(fc|.) 



(14) 



(the coefficients rms and t^s describe scattering into 
evanescent modes). To leading order in t/tc, we recognize 
the standard expressions for the one-dimensional bar- 
rier problem where the reflection and transmission am- 
plitudes can be expressed in terms of momentum ratios. 
This result comes about since, for t = 0, the two phases, 
Sp and the Mott insulator, are essentially identical. How- 
ever, for the current setup, the height of the "barrier" is 
not a free parameter due to the generic nature of the in- 
terface; the result then is determined by the dispersions 
at the boundary. For a small but finite hopping t/tc, the 
Mott insulator remains essentially unchanged, while the 
bosons in the superfluid exploit the kinetic energy. This 
leads to a reduced transmission due to the wave func- 
tion mismatch at the boundary, as reflected by the term 
oc t/tc in Tps. 

Next, we discuss the situation near the tip of the Mott 
lobe; in order to formulate the results for t —^ tc, we first 
define the reduced distance dt = 1 — t/tc > from the 
tip of the Mott lobe and the coherence functions 



f±{k) = [A{k)±iB{k)]i 



t = ta 



(15) 



While the phases arg[/±(fc)] run from arg[/±(0)] = T7''/4 
to arg[/±(7r)] = ±acos(l/\/3 + 1/VG), both moduli di- 
verge at fc — > 0, reflecting the softening of all modes at the 
tip of the Mott lobe; however, the ratio |/+(fc)//_(fc)p = 
1 remains finite. Away from fc — *■ the moduli of f±{k) 
are well behaved. To loading order in 5t we obtain 



.W 



ki 



V2kl 



Ths (e) ~ - 



p ' 



f-jkl) 



f-iM) 



Uiki) f^{ki) 



0{5t), (16) 
'0{5t). 



For 5t 



0, we have kl 



k'l, the coherence fac- 



tors are approximately unity, and we again recognize the 
typical result for a simple barrier. This time, however, 
the particle and the hole branches are degenerate and 
are transmitted equally to leading order, explaining the 
factor \/2. Furthermore, the massive mode is orthogonal 
to the sound mode and does not contribute in leading 
order. Going to nonzero St, the superfluid develops "par- 
ticle" nature, leading to an increase (reduction) in Tps 

(Ths)- 

The behavior of the scattering coefficients for inter- 
mediate values of t/tc is shown in Fig. [2l The coeffi- 
cients are given for an energy e — fkUsiT^)/^ allowing for 
a comparison at different values of t/tc- Note that for 
this energy, the hole and massive modes become non- 



propagating for t/tc "^ 0.971. The corresponding density 
of states effect is responsible for the irregular behavior 
around t/tc ~ 0.971. At small t/tc^ the approximations 
[Eq. (HH)] are in good agreement with the exact values; 
the same is true near the tip of the Mott lobe, where 
the coefhcients [Eq. (fT6|) ] agree well with the exact nu- 
merical values. Here, however, the close vicinity of the 
point where the massive mode and the hole turn evanes- 
cent spoils the applicability of the analytical results much 
faster. 

Given the results above, we note the following generic 
trends: Scattering coefficients relating modes of equal 
type, e.g., particle-type sound in the Sp and particle ex- 
citations in the Mott insulator are large, while scatter- 
ing between unequal modes, e.g., particle- type sound in 
the Sp and hole excitations in the Mott insulator is sup- 
pressed. Furthermore, the energy-dependent point on the 
t/tc axis where the massive and hole modes turn evanes- 
cent controls the breakdown of the approximate formulas 
(fT4|) and (fTBjl . The same statements hold for the respec- 
tive cases of an incoming massive mode and the inverted 
situations at the lower boundary of the Mott lobe, i.e., 
with initial states in S^. 



IV. HEAT CONDUCTIVITY 

As an application, we calculate the transport of heat 
through a Mott-insulating layer within a wedding cake 
structure. The derived heat conductivity k provides in- 
sight into the thermalization process when the lattice 
potential is ramped up,^^ and quantifies to what extent 
two superfluid shells are in thermal contact through the 
Mott layer. In order to calculate k, we consider the heat 
current^! 



Q,/3 •' 



(17) 



where the sum is running over input (a) and output (/3) 
channels. In Eq. p7|) . Na{e) is the density of states in the 
input channels a, v^ (e) is the velocity associated with the 
excitations at energy e in the output channels /3, and we 
integrate over all energies e. The Bose-Einstein distribu- 
tion g{e, T) = [e-x^{e / k-B,T) — 1]^^ controls the occupation 
of the bosonic modes and Tq^ are the transmission ampli- 
tudes connecting input and output channels. The linear- 
response expression Eq. p7)) describes the situation close 
to equilibrium. In a real experiment, the initial state 
after ramping of the optical potential may be far away 
from the bosonic equilibrium distribution g(e, T). Never- 
theless, the calculation of the linear-response result [Eq. 
([T7])] provides some generic insights into the behavior of 
the system which applies to such a situation as well. 

It is instructive to calculate the maximal heat conduc- 
tivity of a homogeneous Mott region. For temperatures 
higher than all energy scales in the problem (the repul- 
sion [/), but lower than the band gap to the next Bloch 



band, the heat conductivity is finite and given by 

r'" dk 



..-^ /■'^Z" dk 
C"= E / T77««(fc)^aW^T5(e,T»f/) 

a—p^h ' 

""" c), (18) 



-/cb (VTTWc - Vi - t/tc 



where we explicitly accounted for the lattice constant a. 
From Eq. (|18p we learn that the coefficient k describes 
the transport of entropy "fce" with velocity Ua/h. Below, 
we take 



\t 



tc) = V2--k^ 
h 



(19) 



as our reference value for the heat conductivity. 

We calculate the heat current [Eq. ([T7])] from one su- 
perffuid shell, point A in Fig. UK a) , to point B in the next 
shell and take the derivative with respect to the temper- 
ature T to obtain k{T) (cf. Fig [5]). The input channels a 
are given by the phonons and massive particles impinging 
onto the Mott insulating phase at point A, i.e., from the 
particle- like superfluid Sp. The output channels on the 
other side of the insulating region are the corresponding 
modes in S/j. 

To obtain the transmission amplitudes Tap, we ap- 
ply a transfer matrix formalism^- to divide the problem 
into three parts, the scattering across the interface be- 
tween the superfluid Sp and the Mott region, point A in 
Fig-QJa), the propagation within the Mott layer, and the 
scattering at the second interface connecting the Mott in- 
sulator and the superfluid S/j [cf. point B in Fig. [TJa)]. 
The scattering at the two interfaces is handled as de- 
scribed in Sec. IIIII above; in addition, scattering ampli- 
tudes into and from evanescent modes now have to be 
accounted for. 

To describe the transfer from the superfluid phase into 
the Mott insulator and vice versa, we match the wave 
functions [Eq. (fT51) ] and their derivatives at the boundary 
(placed at 2/ = 0): we define the wave functions 

and impose the matching conditions ij}^^^{Q) — ■»/'s(0) and 
3yi/'„j(0) = dyipg{0) for the two first components; the 
conditions for the other two components then are fulfilled 
automatically due to T-symmetry. The subscripts I and 
r denote right- and left-moving excitations, respectively. 
This procedure provides us with four relations connecting 
the four amplitudes on one side of the interface with the 
four on the other. Solving for the coefficients Pi(r) and 
Huj.\, we obtain the transfer matrix Mg.Mi defined as 



(Pr\ 
Pi 

\HiJ 
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FIG. 3: Sketch of a scattering event of an incoming phonon from Sp. (a) The incident sound mode is transformed into a 
particle and a hole with amplitudes Tsp and Tsh, respectively. Inside the Mott layer, the particle propagates until it hits the 
potential, where it is either reflected or tunnels under the barrier. The hole emerges from under the barrier and propagates 
to the other end of the layer. At the boundary to Sh, both modes are converted into the modes of the superfluid with the 
corresponding amplitudes, (b) The same situation with a higher energy of the incoming sound mode that allows for undamped 
propagation through the Mott region. 



This transfer matrix depends on the type of interface, Sp 
or S/i, connecting to the Mott layer; correspondingly, we 
denote the two different matrices by MgL'^Mi- Note that 
the four-dimensional character of Eq. (I20|) is due to the 
restriction to the "+" -sector of T; in general, one expects 
the transfer matrix to act on a vector space of twice the 
dimension of the spinor. However, the matrix elements 
connecting the two sectors "-I-" and "— " vanish for a T- 
symmetric system. 

The new task to analyze is the propagation of particle- 
and hole excitations through the Mott region < y < L 
as described by the transfer matrix 



Piiv^L) 

Hriy = L) 

\Hi{v = L)) 



fPr{y^o)\ 
Pi{y = 0) 
Hr{y = 0) 

\Hi{y = Q)j 



(21) 



This task requires the evaluation of the WKB-phases (be- 
tween arbitrary points a and b) 



p{h) 



dxkl^h){x) 



(22) 



into the MI and turn into propagating modes at the right 
of the Mott region. Such turning points are dealt with 
in the standard way-'^ and lead to additional scattering 
phases ±7r/4 in the propagator. 

Depending on the appearance of turning points in the 
particle and/or hole channel, the transfer matrix Mmi as- 
sumes different forms. The simplest case is realized near 
the tip of the Mott lobe, where particle and hole modes 
can propagate unhindered through the Mott region [cf. 
Fig. 2]. In this case, the matrix is diagonal and multi- 
plies each component of the four-spinor with appropriate 
phase factors, 
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The appearance of a turning point in the particle (hole) 
channel renders the 2x2 block matrix m^j (mjjj,) describ- 
ing particle (hole) propagation non-diagonal; note that 
particles (holes) then are reflected but never converted 
into one another. Assuming a classical turning point at 
r^i (cf. trajectory 7^ in Fig. [3]), the transfer matrix for a 
particle excitation assumes the form 



, e 






t3 r n 



-_+i7r/4^ 



of Eq. ID) in the presence of an inhomogeneous poten- 
tial F(y). Attention has to be paied to properly treat 
the classical turning points where the quasi-classical ap- 
proximation [Eq. Il])] breaks down; Fig. [3] illustrates typ- 
ical situations in the present context, where particles and 
holes incident from the left are stopped by the potential 
V^i^ and turn evanescent or tunnel as evanescent modes The equivalent expression for the hole excitation (trajec- 
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FIG. 4: (color online) Left: sketch of the local excitation spectra at different points A, a, fo, c, _B (columns) in the inhomogeneous 
system for three different values of the hopping amplitude at tjtc ~ 0.9, 0.78, 0.25 (rows); Hole- like spectra are dashed; particle- 
like spectra are solid lines. Right: local phase diagram with Mott insulating and superfluid regions depending on the chemical 
potential S^i and the hopping f/fc- The blue (light gray) shaded area of the Mott lobe corresponds to values of t/tc for which 
the hole at point a has spectral overlap with the hole-like spectra at c, allowing for undamped propagation in the energy window 
given in blue (light gray) on the left. The red (dark gray) shaded area of the Mott lobe marks values of t for which the modes 
in the superfluid have mutual spectral overlap for an non-vanishing range of energy as shown on the left. For values t < t, the 
crossed terms do not contribute (phonon ^^ massive mode) to the heat conductivity (see text). 



tory 7ft in Fig. [3]), reads as 



gi(^Q^+i7r/4 



c^'Po,'. 



-i^j. 



-i-KJA 






Note, that due to particle-hole symmetry, the turning 
points r„, are the same for particles and holes. The 
"phases" <^^ {, have to be calculated numerically. The 
total transfer matrix, connecting the two superfluids on 
either side of the Mott insulator, is given by 



M.„ 



(MLmi)" 



(23) 



We obtain the scattering matrix S by solving the linear 
equation 



out 

St 
V Ml'^ J 



MP 



(24) 



for the amplitudes of the outgoing excitations 5^^^, M^^^, 
S^^^, and M^^^. For the heat conductivity, we are only 
interested in the transmission amplitudes connecting the 
channels of incoming phonons and massive modes from 
Sp to the outgoing excitations in S^. Within this sub- 
space the scattering matrix reads as 



(25) 



nis ' mm 






and provides us with the explicit form of the scattering 
amplitudes Tq^ appearing in the expression for the heat 
current [Eq. (fT7)) ]. 

The density of states Na{e) and velocity v 13(e) de- 
rive from the dispersion relations uJs[m) in Eq- ©■ In 
one dimension, these quantities are related via Vs(m) — 
l/-^5(m)(c) &nd hence only in processes where the phonon 
and massive mode channels are mixed with a 7^ /3; there 
appears a ratio Ns(^jn){^) / Nm{s){^) , otherwise all density 
of states effects cancel out. 

Next, we identify those regions in the Mott lobe where 
we expect a large value of k. Figure [5] shows the evo- 
lution of the spectra upon crossing the Mott-insulating 
region, starting with sound- and massive modes in the 
particle-type superfluid Sp at point A right at the inter- 
face, the swap of particle- and hole-type branches within 
the MI regime with a potential rising approximately lin- 
early with distance; see diagrams "a" , "b" , and "c" in 
Fig. |4l and the interchanged massive- and sound modes 
in the hole- type superfluid Sp at point B, again right 
at the interface. Following the evolution of the spectra 
from the tip of the Mott lobe at t = tc down to t = 0, 
various regimes can be identified where transport is fa- 
vored, either via full propagation through the Mott re- 
gion or via conservation of the particle/hole nature of 
the excitation along the trajectory. At large values of 
t close to the tip of the lobe an appreciable part of the 
particle- and hole branches overlaps, allowing these exci- 
tations to propagate through the Mott insulating without 
damping. This overlap terminates when the bottom of 



the hole band lines up with the top of the particle band 
[hu}p{-n/a) — hujh{0)', cf. Eq. ^ and diagram "a" in Fig. 
13], defining the special value i* = (4/5)ic- Another rele- 
vant point is t, where the massive and sound modes stop 
overlapping in the superfluid; see diagram "A" in Fig. |4l 
Comparing the bottom of the massive mode at fc = with 
the top of the sound mode at fc = tt, we find that this 
overlap persists as long as t > i with t = (3 — \/5)ic; 



and we have A,, 



< W^, where Ws 



fiWs{T^/a) de- 



notes the width of the sound mode in Sp. In this sit- 
uation, the particle-type modes in the left (sound) and 
right (massive) superfluids overlap (and vice versa for the 
hole modes) and large transmissions at the boundaries 
enhance the contribution of these modes to k. Finally, 
for t <t, propagation through the Mott region is always 
damped and excitations have to be converted between 
particle- and hole-types, hence only direct terms ex Ir^sp 
(particle- type sound is converted to hole-type sound) and 
oc \Tmm\'^ (hole-typc Higgs is converted to particle- type 
Higgs) contribute. Accordingly, the heat conductivity k, 
naturally splits into direct and cross terms, with the lat- 
ter ones contributing with large weight but only for t > t 
where Ws > A^, 



k{T) = K,„{T) + eiWs - A™)K„_(r), (26) 



with 



Koro 



l^dee\rss\'^+J^dee\r^^f^, (27) 



/A 



For the direct terms, the initial and the final states are 
both sound modes or both massive modes and the density 
of states cancels against the velocity factor. 

The expressions in Eqs. (|27p and ([28]) have to be calcu- 
lated numerically and provide the final result for the heat 
conductivity k{T) shown in Fig. [5] as a function of T/U 
for different values of t/tc- The overall shape involves an 
exponential suppression at small temperatures T, hint- 
ing to the presence of an effective gap A^n, and a sat- 
uration value Koo at large temperatures when all modes 
within the finite bandwidth contribute to the transport. 
Furthermore, the transport efficiency decreases with de- 
creasing t, as is to be expected on the basis of the above 
analysis; cf. Fig. H) 

Given the complexity of the result expressed through 
Eqs. ([?f|) and ([^5]) and the simplicity of the final be- 
havior of K, we have attempted to extract a simple 
and useful expression interpolating between the expo- 
nential rise and the saturation at low and high temper- 
atures. A phenomenological Ansatz with a density of 
states N{e) = S{e — A„h) and a corresponding velocity vq 
provides us with the simple formula 



k(T) w k^voiA^ti) 



/A„„\ exp(A„„/fcBr) 



\kBTj [exp(A„„/fcBr) - 1]2 



(29) 
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FIG. 5: (a) Heat conductivity k{T) for different values of 
t {t = 0.9, 0.82, 0.78,0. 74tc, from the back to the front) for 
a Mott layer of thickness Lq — 10a a,t t — {a denotes the 
lattice constant). The structure of k{T) for small T is dictated 
by the spectral gap in the intermediate Mott region. For 
smaller values of t, where only damped modes contribute to 
transport, the heat conductivity is exponentially suppressed 
and not shown here, (b) The dependence of the maximal heat 
conductivity Ktx> = k{T ^ U/ks) on the hopping amplitude t 
displaying the qualitative change for t < t in the plot, where 
the strong suppression sets in. The change at t* is masked by 
a scattering resonance in the intermediate t <t <ti, regime. 



with x^ exp(a;)/[exp(2:) — Ij^l^^^o ~* 1; we find the ve- 
locity Vq parameter related to the saturation value of k 
at large temperatures, k^o = k^v^; the scaled function 
Koo {t) shown in Fig. [5I^b) reproduces the expected quali- 
tative behavior, with a shoulder at i « t and a suppres- 
sion for t < t^,. The transmittance of the boundaries as 
well as scattering resonances within the Mott layer (be- 
tween the two boundaries to the superfluids or between 
one such boundary and a classical turning point) arc all 
absorbed in the coefficient Kqo- Also, we note that the 
band edges in the excitation spectra manifest themselves 
in the densities of states and give rise to sharp features 
in Koo- These effects are masked in an experiment, where 
the finite size of the superfluid and Mott-insulating re- 
gions induces an uncertainty in the momenta of all exci- 
tations and we account for this smearing in our numerical 
evaluation of Eqs. ^ and (P5| .^^ 

Close to the tip, where the conduction is dominated by 
itinerant modes, the effective gap parameter A,fc turns 
out to be independent of the Mott layer thickness L and 
is approximately given by 



tc-i + 0.32. 



(30) 



The linear behavior in t of A^^ has to be compared to 
the evolution of the size of the maximal gap Am within 
the Mott layer which has a square root dependence on 
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1 — t/tc- The different scaling suggests that the gap A^,- 
plays the role of an effective parameter describing the 
complex transport involving all of the interfaces and the 
inhomogeneous Mott layer and is not directly related to 
the spectral gap in the insulating region. In the same 
way, the offset by ^ 0.32 effectively accounts for the con- 
version of quasi-particles at the phase boundaries. With 
the above choices for the two phenomenological param- 
eters vq and Aef£, we find excellent agreement between 
the numerical data and the results of our simple Ansatz. 
In the regime t < t, the effective gap depends on the 
thickness Lq of the Mott layer (Lq is taken at i = 0). In 
Fig. [5l we show the results for Lq — 10a. The slope of 
the exponential decrease of Koo for t < i depends strongly 
on Lq. However, its significant suppression at Lq = 10a 
indicates that already at this width, the superfluid layers 
are essentially decoupled even at temperatures k^^T w U. 



CONCLUSIONS 



method. 

For a phonon incident from a Sp onto a Mott insula- 
tor, we find standard expressions for the scattering am- 
plitudes describing the scattering at a potential barrier 
in the limits i — *■ 0, tc, with the involved momenta deter- 
mined by the non-trivial bulk dispersions. Going away 
from the Mott-lobe base and tip, the amplitudes pick 
up non-trivial corrections which are easily found numer- 
ically. 

In calculating the heat conductivity across a Mott layer 
in a wedding cake structure, we have combined the scat- 
tering at the two interfaces with the quasi-classical prop- 
agation through the inhomogeneous Mott layer. We find 
that for a Mott shell at moderately to large hopping, i.e., 
t > O.Stc, the adjacent superfluid shells are in good ther- 
mal contact. For small hopping below t, however, the 
Mott shells represent practically infinite barriers. Impli- 
cations of our findings on the lattice ramping probleroii 
deserve further studies. 



Summarizing, we have developed a framework to ad- 
dress the behavior of quasi-particle excitations in a 
strongly correlated bosonic heterostructure. We have 
derived a set of first-quantized spinor wave functions 
valid in the superfluid- and Mott-insulating phases, have 
derived the scattering properties of a superfluid-Mott- 
insulator interface, and have calculated the heat con- 
ductance across a Mott region as an application of our 



Acknowledgments 

We thank F. Hassler and L. PoUet for extensive dis- 
cussions and acknowledge financial support from the 
Swiss National Foundation through the NCCR MaNEP. 
S.D.H. acknowledges the hospitality of the Institute 
Henri Poincare-Centre Emile Borel. 



D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and 
P. Zoller, Cold bosonic atoms in optical lattices, Phys. Rev. 
Lett. 81, 3108 (1998), URL, 

^ M. Lewenstein, A. Sanpera, V. Ahufinger, D. B., 
A. Sen De, and U. Sen, Ultracold atomic gases in optical 
lattices: mimicking condensed matter physics and beyond, 
Adv. in Phys. 56, 243 (2007), URL 

^ S. Wessel, F. Alet, M. Troyer, and G. G. Batrouni, Quan- 
tum Monte Carlo simulations of confined bosonic atoms in 
optical lattices, Phys. Rev. A 70, 053615 (2004), URL 

" O. Gygi, H. G. Katzgraber, M. Troyer, S. Wessel, and G. G. 
Batrouni, Simulations of ultracold bosonic atoms in optical 
lattices with anharmonic traps, Phys. Rev. A 73, 063606 

(2006), URL, 

® V. A. Kashurnikov, N. V. Prokof'ev, and B. V. Svistunov, 
Revealing the superfluid-Mott-msulator transition in an 
optical lattice, Phys. Rev. A 66, 031601(R) (2002),!^^ 

^ S. Foiling, A. Widera, T. Miiller, F. Gerbier, and I. Bloch, 
Formation of Spatial Shell Structure in the Superfluid to 
Mott Insulator Transition, Phys. Rev. Lett. 97, 060403 
(2006), URL 

"^ M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. 
Fisher, Boson localization and the superfluid-insulator 
transition, Phys. Rev. B 40, 546 (1989), URL 

* A. F. Andreev, The Thermal Conductivity of the Interme- 
diate State in Superconductors, Sov. Phys. JETP 19, 1228 
(1964). 



® S. Vishveshwara and C. Lannert, Josephson physics medi- 
ated by the Mott insulating phase, Phys. Rev. A 78, 053620 
(2008), ^D 

^° S. D. Huber, E. Altman, H. P. Biichler, and G. Blatter, Dy- 
namical properties of ultracold bosons in an optical lattice, 
Phys. Rev. B 75, 085106 (2007), URL, 

" L. PoUet, C. KoUath, K. van Houcke, and M. Troyer, Tem- 
perature changes when adiabatically ramping up an optical 
lattice. New J. Phys. 10, 065001 (2008), URL, 

^^ C. Honerkamp and M. Sigrist, Andreev Reflection in Uni- 
tary and Non-Unitary Triplet States, J. Low Temp. Phys. 
Ill, 895 (1998), URL 

^^ G. E. Blonder, M. Tinkham, and T. M. Klapwijk, Transi- 
tion from metallic to tunneling regimes in superconducting 
microconstrictions: Excess current, charge imbalance, and 
supercurrent conversion, Phys. Rev. B 25, 4515 (1982), 
URL 

^'^ A. B. Migdal, Qualitative Methods in Quantum Theory 
(Advanced Books Classics, 1977). 

^^ G. Wentzel, Eine Verallgememerung der Quantenbedin- 
gungen fiir die Zwecke der Wellenmechanik, Z. Phys. A 
38, 518 (1926), URL 

^^ H. A. Kramers, Wellenmechanik und halbzahlige Quan- 
tisierung, Z. Phys. A 39, 828 (1926), URL, 
L. Brillouin, La mecanique ondulatoire de Schrodinger; une 
methode generale de resolution par approximations succes- 
sives, C. R. Acad. Sci. 183, 24 (1926). 



11 



M. Y. Azbel, Eigenstates and properties of random systems 

in one dimension at zero temperature, Phys. Rev. B 28, 

4106 (1983), lURL] 

This would not be the case for a discussion of the Josephson 

effect. 

C. Kollath, U. Schollwock, and W. Zwerger, Spin-Charge 

Separation in Cold Fermi Gases: A Real Time Analysis, 

Phys. Rev. Lett. 95, 176401 (2005), tJRLl 



G. D. Mahan, Many-Particle Physics (Plenium Press, New 
York and London, 1990), p. 31. 

In order to account for this smearing, we cut off the in- 
tegrals at low momenta kc ex. 1/Lo, and add a convo- 
lution to the integrals in Eq. (|27|l and (|28p of the form 
J de ^ J dede'f{e — e') where f{x) is a Lorentzian which 
is centered at a; = and has a width oc 1/Lq. 



